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Abstract 

In general the calculation of robustness of entanglement for the mixed entangled 
quantum states is rather difficult to handle analytically. Using the the convex semi- 
definite programming method, the robustness of entanglement of some mixed entangled 
quantum states such as: 2 ® 2 Bell decomposable (BD) states, a generic two qubit state 
in Wootters basis, iso-concurrence decomposable states, 2 3 Bell decomposable states, 
d® d Werner and isotropic states, a one parameter 3 ® 3 state and finally multi partite 
isotropic state, is calculated exactly, where thus obtained results are in agreement with 
those of :2 (g> 2 density matrices, already calculated by one of the authors in [1, 2]. Also 
an analytic expression is given for separable states that wipe out all entanglement and 
it is further shown that they are on the boundary of separable states as pointed out in 
[3]. Keywords: Robustness of entanglement, Semi-definite programming, Bell 
decomposable states, Werner and isotropic states. 

PACs Index: 03.65.Ud 
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Quantum entanglement has recently been attracted much attention as a potential resource for 
communication and information processing [4, 5]. Entanglement usually arises from quantum 
correlations between separated subsystems which can not be created by local actions on each 
subsystem. By definition, a mixed state p of a bipartite system is said to be separable (non 
entangled) if it can be written as a convex combination of product states 

i i 

where pf^ and pf^ are states of subsystems 1 and 2, respectively. Although, in the case of 
pure states of bipartite systems it is easy to check whether a given state is, or is not entangled, 
the question is yet an open problem in the case of mixed states. Therefore having a measure 
to quantify entanglement of mixed states is likely to be valuable and several measures of 
entanglement have been proposed [6, 7, 8, 9]. 

One useful quantity introduced in [10] as a measure of entanglement is robustness of entan- 
glement. It corresponds to the minimal amount of mixing with locally prepared states which 
washes out all entanglement. An analytical expression for pure states of bipartite systems by 
using partial transpose has been given in [10] . Meanwhile the robustness of entanglement has 
been calculated for a Werner states. Moreover, in [11] Vidal and Werner have computed the 
robustness of entanglement for density operators with symmetry. In [12] Rudolph, using cross 
norm has clarified the relationship of the greatest cross norm with the robustness of entan- 
glement and has determined the value of the greatest cross norm for Bell diagonal states. A 
geometrical interpretation of robustness is given in [3] and it is pointed out that two corre- 
sponding separable states needed to wipe out all entanglement are necessarily on the boundary 
of separable set. On the other hand, the robustness of entanglement of few mixed quantum 
states such as: 2 ® 2 Bell decomposable (BD) states and a generic two qubit state in Wootters 
basis is already calculated by one of the authors in[l, 2]. In Ref [13] has characterized the 
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robustness of entanglement, and its relation to the permutation symmetries, for the basic set 
of eight entangled three particle states of spin-1/2 objects. Authors in [14], have studied the 
robustness of multi-party entanglement under local decoherence, modeled by partially depolar- 
izing channels acting independently on each subsystem. Unfortunately, in general, the above 
mentioned quantity as the most proposed measures of entanglement involves exteremization 
which is difficult to handel analytically. 

On the other hand, over the past years, semidefinite programming (SDP) has been rec- 
ognized as valuable numerical tools for control system analysis and design. In (SDP) one 
minimizes a linear function subject to the constraint that an affine combination of symmet- 
ric matrices is positive semidefinite. SDP, has been studied (under various names) as far 
back as the 1940s. Subsequent research in semidefinite programming during the 1990s was 
driven by applications in combinatorial optimization[15], communications and signal process- 
ing [16, 17, 18], and other areas of engineering[19]. Although semidefinite programming is 
designed to be applied in numerical methods it can be used for analytic computations, too. 
Some authors try to use the SDP to construct an explicit entanglement witness [20, 21]. Ki- 
taev used semidefinite programming duality to prove the impossibility of quantum coin flipping 
[22], and Rains gave bounds on distillable entanglement using semidefinite programming [23]. 
In the context of quantum computation, Barnum, Saks and Szegedy reformulated quantum 
query complexity in terms of a semidefinite program [24]. The problem of finding the opti- 
mal measurement to distinguish between a set of quantum states was first formulated as a 
semidefinite program in 1972 by Holevo4, who gave optimality conditions equivalent to the 
complementary slackness conditions [25]. Recently, Eldar, Megretski and Verghese showed 
that the optimal measurements can be found efficiently by solving the dual followed by the 
use of linear programming [26] . Also in [27] used semidefinite programming to show that the 
standard algorithm implements the optimal set of measurements. All of the above mentioned 
applications indicate that the method of SDP is very useful. 
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Here in this paper, by using the the convex semi-definite programming method, the robust- 
ness of entanglement of some mixed entangled quantum states such as: 2(g) 2 Bell decomposable 
(BD) states, a generic two qubit state in Wootters basis, iso-concurrence decomposable states, 
2 ® 3 Bell decomposable states, d <S> d Werner and isotropic states, a one parameter 3 <8> 3 state 
and finally multi partite isotropic state, is calculated exactly, where thus obtained results are 
in agreement with those of: 2 ® 2 density matrices, already calculated by one of the authors 
in[l, 2]. Also an analytic expression is given for separable states that wipe out all entanglement 
and it is further shown that they are on the boundary of separable states as pointed out in [3] . 

The paper is organized as follows: 
In sections 2 and 3 we give brief review of semidefinite programming and robustness of entan- 
glement, respectively. In section 4, by using the semi-definite programing method we calculate 
the robustness of entanglement of some mixed entangled quantum states, such as: 2 <g> 2 Bell 
decomposable (BD) states, a generic two qubit state in Wootters basis, iso-concurrence de- 
composable states, 2 <g> 3 Bell decomposable states, d <g> d Werner and isotropic states, a one 
parameter 3 ® 3 state and finally multi partite isotropic state. The paper is ended with a brief 
conclusion. 

2 Semi-definite programming 

A semidefinite programming(SDP) is a particular type of convex optimization problem [28]. 
A semidefinite programming problem requires minimizing a linear function subject to a linear 
matrix inequality (LMI) constraint [29]: 

minimize V = c T x 

(2-1) 

subject to F(x) > 0, 

where c is a given vector, x = (x±, ...,x n ), and F(x) = F + J2i%iF, for some fixed hermitian 
matrices Fj. The inequality sign in F(x) > means that F(x) is positive semidefinite. 
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This problem is called the primal problem. Vectors x that satisfy the constraint F(x) < 
are called primal feasible points, and if they satisfy F(x) > they are called strictly feasible 
points. The minimal objective value c T x is by convention denoted as V* and is called the 
primal optimal value. 

The minimization is performed over the vector x, whose component are the variables of the 
problem. The vector x which satisfies the LMI, is called a feasible solution, and the set of all 
feasible solutions, is called the feasible set. 

A very important property of a (SDP) is its convexity, since the feasible set defined by 
the above constraints is convex. For this reason, semidefinite programming has a nice duality 
structure, with, the associated dual program being: 

maximize —Tr [F Z] 

Z > (2-2) 
Tr[FZ] = a. 

Here the variable is the real symmetric (or Hermitean) matrix Z, and the data c, Fi are 
the same as in the primal problem. Correspondingly, matrices Z satisfying the constraints are 
called dual feasible (or strictly dual feasible if Z > 0). The maximal objective value —TrF Z, 
the dual optimal value, is denoted as d*. 

The objective value of a primal feasible point is an upper bound on V*, and the objective 
value of a dual feasible point is a lower bound on d*. The main reason why one is interested in 
the dual problem is that one can prove that, under relatively mild assumptions, V* = d*. This 
holds, for example, if either the primal problem or the dual problem are strictly feasible, i.e. 
there either exist strictly primal feasible points or strictly dual feasible points. If this or other 
conditions are not fulfilled, we still have that d* < V*. Furthermore, when both the primal 
and dual problem are strictly feasible, one proves the following optimality condition on x: 

A primal feasible x and a dual feasible Z are optimal which is denoted by x and Z if and 
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only if 

F{x)Z = ZF{x) = 0. (2-3) 

This latter condition is called the complementary slackness condition. 

In one way or another, numerical methods for solving SDP problems always exploit the 
inequality d < d* < V* < V, where d and V are the objective values for any dual feasible 
point and primal feasible point, respectively. The difference 

V - d = c T x + Tr[F Z] = Tr[F(x)Z] > (2-4) 

is called the duality gap, and the optimal value V* is always bracketed inside the interval 
[d, V] . These numerical methods try to minimize the duality gap by subsequently choosing 
better feasible points. Under the requirements of the above-mentioned theorem, the duality 
gap can be made arbitrarily small (as far as numerical precision allows). 

Equation (2-3) together with (2-2) and (2-1) constitute a set of necessary and sufficient 
conditions for x to be an optimal solution to the problem of (2-1), when both the primal and 
the dual are strictly feasible. 

3 Robustness of entanglement 

According to [10] for a given entangled state p and separable state p", a new density matrix 
p' s ) can be constructed as, 

p' s = -^—(p + sp"), s>0, (3-5) 
s + 1 

where it can be either entangled or separable. It was pointed that there always exits the 
minimal s corresponding to one p" such that p' s is separable. This minimal s is called the 
robustness of p relative to p", denoted by R(p\\p"). The absolute robustness of p is defined as 
the quantity, 

R(p\\S) = minR(p\\p"), p" G S, (3-6) 
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where S is the set of separable states. 

Du et al. in [3] have given a geometrical interpretation of robustness and pointed that if 
s in Eq. (3-5) is minimal among all separable states p" , i.e. s is the absolute robustness of p, 
then p" and p' s in Eq. (3-5) are necessarily on the boundary of the separable states. 

4 Robustness of entanglement via semi-definite program- 
ming 

Unfortunately, the above mentioned quantity as the most proposed measure of entanglement 
involves extremization which are difficult to handle analytically. One of authors have given 
analytical expression for the robustness of entanglement of some 2® 2 density matrices in [1, 2], 
here in this section we try to obtain robustness of entanglement for many categories of states, 
namely, 2 ® 2 Bell decomposable (BD) states, a generic two qubit state in Wootters basis, iso- 
concurrence decomposable states, 2 <g> 3 Bell decomposable states, d <S> d Werner and isotropic 
states, a one parameter 3 <8> 3 state and finally multi partite isotropic state, via semi-definite 
programming method. As we will see in the following, besides the elegance of semi-definite 
programming in the calculation of the robustness of entanglement, also we do not need to 
define any kind of norm for mixed quantum states in order to calculate their robustness of 
entanglement as it is done in Ref. [1, 2]. 

4.1 Robustness of entanglement for Bell-decomposable state 

A Bell decomposable (BD) state is defined by: 



4 



4 



(4-7) 



i=i 

where is Bell state, given by: 




(4-8) 
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M = 



V2 



(ITT) - III)), 



m= ^ + ) = ^ ! (iti) + iit)), 



lifc> = iT) = ^(IU> - UT»- 



In terms of Pauli's matrices, p can be written as, 



1 3 

p = -(I <g) / + ^tiCTj <g) CTj), 
4 i=l 



where 



(4-9) 
(4-10) 
(4-11) 



(4-12) 



From the positivity of p we get 



h = Pi -P2+P3-P4, 

*3 = Pl+P2-P3~Pi- 
1 + tl - f 2 + t 3 > 0, 

1 - tl + t 2 + h > 0, 

1 + h + 1 2 - 1 3 > 0, 

1 - tx - t 2 - t 3 > 0. 



(4-13) 



(4-14) 



These equations form a tetrahedral with its vertices located at (1, —1, 1), (—1, 1, 1), (1, 1, —1), 
(—1,-1,-1) [30]. In fact these vertices denote the Bell states given in Eqs. (4-8) to (4-11), 
respectively. 

On the other hand p given in Eq. (4-12) is separable if and only if U satisfy Eq. (4-14) 
and, 

1 + h + 1 2 + 1 3 > 0, 

1 - *1 - *2 + *3 > 0, 
1 + f 1 - t 2 - *3 > 0, 

1 - h + 1 2 - u > 0. 



(4-15) 
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Inequalities (4-14) and (4-15) form an octahedral with its vertices located at Of = (±1, 0, 0), 
Of = (0, ±1,0) and Of = (0,0, ±1). So, tetrahedral is divided into five regions. Central 
regions, defined by octahedral, are separable states (p k < |). There are also four smaller 
equivalent tetrahedral corresponding to entangled states(pfc > \ for only one of k — 1, ...,4), 
where Pk = \ denote to boundary between separable and entangled region. Each tetrahedral 
takes one Bell state as one of its vertices. Three other vertices of each tetrahedral form a 
triangle which is its common face with octahedral (See Fig. 1). 

Here in this section we evaluate robustness of entanglement for all BD-states with semi- 
definite programming method, and we give an explicit form the corresponding p' s and p" which 
are on the boundary of the separable states. 

Now for a given BD density matrix 

4 1 4 1 

i=l 1 i=2 1 

and arbitrary separable density matrix 

4 i 

p1 = £p1|^><^|, with P ;<-, (4-17) 

1=1 

according to the SDP method explained in section (2), we have to optimize V = c T x = 
—Tr(Ap) with 

F(x) = F + AFi = p' s + ^—{-P) > 0, (4-18) 

1 + s 



Therefore, we have 



V = -A, and A = x = (4-19) 

1 + s 



Now using the complementary slackness equation (2-3) with a optimal feasible pair (Z, A), we 
have 

Z(p' s -kp) = 0, (4-20) 



or 



Z(I-App i s )=0, (4-21) 
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where p' s is pseudo inverse of p' s . By substituting (4-16) and (4-17) into (4-21) and considering 
the positivity of p' s — Ap and after some elementary algebra we arrive at the following results, 



iiii 



A = m m{^A-,-}- (4-22) 
Pi P2 P3 Pa 

Now defining 

sff' = {l + s)p a -p, (4-23) 

according to equation (3-5) we get 

sp'l+Pi = (l + s)p i , i = 1,2,3,4 (4-24) 
Hence using the above equation we get the following result for the parameter 

, = ^ (4-25) 

Pi- P i 

The choice A = ^ leads to the rank three density matrix p" with p'[ = and the parameter 

* = ^ (4-26) 

Pi 

l l 

but other choices of A = P - L , % — 2, 3, 4 leads to rank three matrices with p'l ^ and ss = p '~ p }, . 
Comparing s £ , i = 1,2,3,4, one can show that S\ is smaller than the others, hence for given 
separable density matrix p s , the choice A = ^ yields the minimum parameter s±. 

We see that the parameter si = — 1 is a monotonic decreasing function of p\ in the 
separable region < p\ < 1/2 and its minimum value can be obtained for p[ = 1/2 states 
which lies at the boundary of separable region. Therefore the robustness of entanglement, that 
is, the minimum of s 

s = 2 Pl - 1 (4-27) 

corresponds to separable states p and p" lying at the corresponding boundaries p[ = | and 
p'[ = of separable region, in agreement with the results of references [31, 1, 3]. 

So far using the SDP optimization method we have proved that for a given entangled 
density matrix p minimum s in formula (3-5) is achived for separable states p' s and p" lying 
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at the boundary of separable region. One we choose p" at that part of boundary of separable 
region far from p, one can determine p' s simply from the intersection of a straight line drawn 
from p to p" and the segment of the boundary of separable region near to p. As it is shown 
in Fig (2) the boundary S[ = P1P2P3 is near p and others S 2 = PiP^O^O^ , S 3 = PiP 2 Oi O2 , 
S 4 = P2P3O1O3 are far from it. Therefore, first we have to choose p" at one of S' 1: <S 2 , <S 3 , <5 4 
, but we see that only the choice of S[ leads to p'{ = 0, that is minimum s\ and other choices 
leads to greater value of parameter s. But the choice of p" at any point of ^=boundary will 
yields the same minimum robustness of entanglement, which result in the choice of p' s belong 
to a triangle defined by 

^ - Wl — 2> P2 — 2 4pi ' P3—2 4pi ' P4— 2pi 2/' 

d _ r ' _ 1 ' _ 1 _ i-2p2 ' _ 1-P2-P4 _ 1 ' _ 1 _ 1-2P4 1 (4_28) 

D — \P\ — 2' ^2 - 2 4pi ' P3 — 2 P1 2' ^4 - 2 4 P1 J ' V ; 

/-»_/'_! ' _ 1-P3~P4 _ 1 ' _ 1 _ l-2p 3 ' _ 1 _ l~2p4 1 

^ — \P\ — 2 i Pi— 2 pi 2' ^3— 2 4pi ' ^4— 2 4pi J' 

lying at boundary S[. 

4.2 Robustness of entanglement for 2x2 density matrix in Woot- 
ters's basis 

Here, we find robustness of a generic two qubit density matrix. To this aim we first review 
Wootters's basis as presented by Wootters in [9]. Wootters in [9, ?] has shown that for any 
two qubit density matrix p there always exist a decomposition 

P = £NN, (4-29) 

i 

called Wootters's basis, such that 

(xi\xj) = XiSij, (4-30) 

where Aj are square roots of eigenvalues, in decreasing order, of the non-Hermitian matrix pp 
and 

P = (<Ty <S> Oy)p*((Jy <g) (Ty), (4"3l) 
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where p* is the complex conjugate of p when it is expressed in a standard basis such as 
{ITT) j |t!)}j {|!t) ) I!!)} an d c y represent Pauli matrix in local basis {|T) , I!)} • Based on this, 
the concurrence of the mixed state p is defined by max(0, Ai — A 2 — A 3 — A 4 ) [9] to see the 
explicit form of the wootters basis of the generic 2*2 density matrix see ref . [33] . 

Now let us define states |x-) as 

\x'i) = -^=, for i = 1,2,3,4. (4-32) 



Then p can be expanded as 

P = E A *K)(^I, (4-33) 

i 

and Eq. (4-30) takes the following form 

Vj^) = Sij. (4-34) 

Here in this section we obtain the robustness for a generic two qubit density matrix with 
SDP method. Our method of evaluation of robustness is based on the decomposition of 
density matrix given by Wootters in [9]. By defining Pj = A^-fQ, where K t = (x'^x'j), then 
normalization condition of p leads to 

4 

Tr(p)=]Tp = l, P>0. (4-35) 
i=i 

This means that with respect to coordinates Pj, the space of density matrices forms a tetra- 
hedral. With respect to this representation separability condition Ai — A 2 — A3 — A4 < takes 
the following form 

Pl P * Ps Pa < 0. (4-36) 



K x K 2 K 3 K A 

The states that saturate inequality (4-36) form a plane called «Si (see Fig. 3). All states 
violating inequality (4-36) are entangled states for which Ai > A 2 + A3 + A4. These states 
form an entangled region with <Si as its separable boundary. There exist, however, three other 
entangled regions corresponding to the dominating Xj (j = 2, 3, 4), respectively. These regions 
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also define separable planes Sj. Four planes Si together with four planes <S ? ', corresponding to 
Aj = 0, form an irregular octahedral corresponding to the separable states. This geometry is 
similar to that of Bell decomposable states but here we have an irregular octahedral associated 
to separable states [32]. Figure 3 shows a perspective of this geometry, where two separable 
planes Si and S[ are shown explicitly. 

Now in order to obtain the robustness of p, suppose that a ray from p is drawn such that 
intersects the boundary planes of separable region at points p' s and p", respectively. Although 
p' s is necessarily on the plane S 1 , but p" is allowed to lie on any plane S[, S 2 , S 3 or <S 4 , where 
we evaluate robustness for each case separately. 

Again for a given generic entangled density matrix p in wootters basis 

P = Y, x i\x'i)(x'i\, Xi > A 2 + A 3 + A 4 (4-37) 
and an arbitrary separable density matrix in the same wootters basis 

p s = J2 x 'i\ x 'i)( x 'ii a;<a; + a; + a; (4-38) 

the SDP optimization of — Tr(Ap) with respect to p s — Ap > yields 



A = min{± ±, ^ ^}. (4-39) 

M ^2 ^3 ^4 

In this case p" can be written as a convex sum of three vertices of the plane 

p" = X i \ x 'i) (41 — °2°"2 + a 3 <T3 + a 4 a 4 , a 2 + a 3 + a 4 = 1, (4-40) 

i 

where <7j are separable states that can be written as a convex sum of two corresponding vertices 
of tetrahedral as 

a 2 = — — - |4) (4| + — -f — |4) (4| , (4-41) 





1 


14) (41 


+ 




1 


14) 


(41 




+ A^ 4 




K 3 


+ A^ 4 




1 


\ X 2) ( X 2\ 


+ 




1 


14) 


(41 


K 2 


+ A^ 4 


K 2 


+ K 4 




1 


\ X 2) ( X 2\ 


+ 




1 


14) 


(41 


K 2 


+ K 3 


K 2 


+ K 3 



a 3 = |4) (4| + |4) (41 , (4-42) 

A 2 + A 4 K 2 + A 4 

a 4 = — -1— |4) (4| + — -1— |4) (4| , (4-43) 
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and A" are 

A" = 0, (4-44) 

X-T^ + T^'xTk^-X"- (4 " 47) 

By expanding p' s as convex sum of p and p" and also using the fact that the coordinates of p' s 
satisfy the equation 

pi pi pi pi 

k-k-k-^r ' (4 - 48) 

after some algebra, coordinates A^ of p' s can be written as 

I ai I a-A I Q4 A \ 

w _ 1,^3+^4 ^ K 2 +K 4 K 2 +K : J A1 , , Hf o, 

a i - a i a i "4 |~' t 4 " 4 ^ 

^3+^4 ^ K2+K4 ^ K 2 +K 3 ^ 2 

( a 2 I Q3 I "4 A \ I C / Q3 I "4 A 

y _ Vi^3+^4 3 X2+K4 3 K2+K3) 7X2 T 2 VX 2 +^4 3 K 2 +K 3 ) (4-501 

2 (12 I 03 1 a 4 1 C ' \ / 

i^3+^4 X 2 +i^4 ^ K 2 +K 3 

( "2 I 03 I __£4__\ \ I C ( 0,2 1 Q4 \ 

w = V^ 3 +^4 _ K 2 +K 4 3 K 2 +K 3 J A 3 T 2 Vx 3 +i^ 4 K 2 +K 3 ) (A-hl) 

3 J22 |_ a 3 I a 4 1 C ' \ / 



J <£3 L 

if 3 +^4 ^ X 2 +i^4 ^ ^2+^3 



Q2 J 03 



J "3 1_ 





+ 


u-4 
K2 + K3 


T 2 






<J2 


+ 




^3+^4 


K 2 +K 3 


T 2 






<J2 


+ 


0.3 



y _ \-g3+^"4 X2+-ftT4 i^2+-^3/ ~* ' ^ yA 3 +A 4 ' A2+A4/ ('4-52") 

4 a2 I a 3 1 a 4 1 C ' \ / 

1^3+^4 ~T X 2 +i^ 4 i^2+X 3 2 

where C = Ai — A 2 — A 3 — A 4 is the concurrence of p. By using the above result one can evaluate 
robustness of p relative to p" as 

1 C 

Sl = T _ 1 = 2a 2 1 2a 3 1 2^~ " (4-53) 

A K 3 +K 4 K 2 +K 4 K 2 +K 3 



Analogue to the above method one can evaluate robustness of p for the case that p" lies on 

A 2 ' 



the plane S 2 i.e. A = ^. In this case p" s can be expanded as convex sum of three vertices of 



the plane 

p" = Mi + + M4, 61 + 63 + 64 = 1, (4-54) 
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where 

o x = 1 \x\) (x\\ + 1 |4) (4| , (4-55) 

and <j 3 and 04 are defined in Eqs. (4-42) and (4-43). Then after some algebra we obtain the 
corresponding robustness as 

C 

S 2 = 2fc3 2b 4 ( 4 - 56 ) 
K2+K4 ^ K 2 +K 3 

Similarly in cases that separable state p" s are on the planes S 3 and S4 we obtain relative 
robustness of p as 

*3 = 2C9 , ° 2c , , (4-57) 

K 3 +K4, ^ K 2 +K 4 

and 

C 

S 4 = 2ds 1 2d, ~ » ( 4 " 58 ) 

respectively. Equations. (4-53), (4-56), (4-57) and (4-58) show that in order to achieve the 
minimum robustness it is enough to consider the case that separable state p" lies on the 
plane S[. Hence the choice A = ^ corresponds to smaller value of parameter s\ than others. 
Therefore for the separable density matrix p' s given in (4-38), the minimum value of parameter 
si, is given by 

ai = ^-l = ^-l (4-59) 
Ai Pi 

Obviously Si is a decreasing function of parameter p' v Again the minimum value of si can be 
obtain from maximum possible value p[, since si is a monotonic decreasing function of p\, and 
its maximum value corresponds to the separable states with 

a'i = a; + a; + a;, (4-eo) 

i.e. , the separable state p' s lying on it boundary of separable region in agreement with [31, 1, 3]. 
With this consideration we are now allowed to choose coefficients in such a way that Eq. 
(4-53) becomes minimum. It is easy to see that this happens as long as the coefficient a& 
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corresponding to the term min(Ki + Kj) becomes one. Therefore the robustness of p relative 
to p" s is 

s = min(i ^ +i ^C (4-61) 

which is one of the main results of this work. Here the minimum is taken over all combination 
of Ki + Kj for i,j = 2,3,4. Equation (4-61) implies that for two qubit systems robustness 
is proportional to the concurrence. We see that the minimum robustness given in Eq. (2-3) 
corresponds to Oj = 5 ik , therefore, by using Eq. (4-40) we get the following result for p" 

Ps = °h- (4-62) 

Also by using = in Eqs. (4-49) to (4-52) one can easily obtain the coordinates A- of 
separable state p' s . 

As we will show in the next section the Bell decomposable states correspond to the Ki = 1 
for i = 1,2,3,4, therefore in Bell decomposable states Eq. (4-61) implies that the robustness 
is equal to the concurrence. 

One can show that thus obtained robustness is minimum over all separable states. 

4.3 Iso-concurrence decomposable states 

In this section we define iso-concurrence decomposable (ICD) states, then we give their sep- 
arability condition and evaluate robustness of entanglement. The iso-concurrence states are 
defined by 

=cos0|OO) +sin0|ll), \4> 2 ) = sin0|OO) -cos0|ll), (4-63) 
|0 3 > = cos0|Ol) + sin0|lO), \<f> 4 ) =sin0|Ol) -cos0|lO). (4-64) 

It is quite easy to see that the above states are orthogonal, thus span the Hilbert space of 2 <g> 2 
systems. Also by choosing = f the above states reduce to Bell states. Now we can define 
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ICD states as 

4 4 

P = £Pi|&X&l, 0< Pi <l, 5> = 1. (4-65) 

i=l i=l 

These states form a four simplex (tetrahedral) with its vertices defined by p\ = 1, p 2 = 1, 
p 3 = 1 and p 4 = 1, respectively. 

Peres-Horodeckis criterion [34, 35] for separability implies that the state given in Eq. (4-65) 
is separable if and only if the following inequalities are satisfied 



(pi - 


-P2) 


< \/4p3P4/sin 2 2^+ (pa 


-P4) 2 , 


(4-66) 


(P2 " 


-Pi) 


< 7^4/ sin 2 2^+ (pa 


-p 4 ) 2 , 


(4-67) 


(P3 




) < v ^4pip 2 sin 2 26 l + (pi 


-P2) 2 , 


(4-68) 


(P4 " 


-Ps) 


< v / 4piP 2 /sin 2 2^+ (px 


-P2) 2 . 


(4-69) 



Inequalities (4-66) to (4-69) divide tetrahedral of density matrices to five regions. Central re- 
gions, defined by the above inequalities, form a deformed octahedral and are separable states. 
In four other regions one of the above inequality will not hold, therefore they represent entan- 
gled states. Bellow we consider entangled states corresponding to the violation of inequality 
(4-66) i.e. the states which satisfy the following inequality 

(Pi - P2) > \/ 4 P3P4/sin 2 29 + (p 3 -p 4 ) 2 . (4-70) 

In order to obtain the robustness of ICD states, we have to follow the method presented 
by Wootters in [9]. Starting from the spectral decomposition for ICD states, and defining 
subnormalized orthogonal eigenvectors, the wootters basis of ICD states can be defined as 



\ x l) = 


—iOL\ 


sfp~\\4>^ + ioc 2 ^ 


/P2I02), 




\x 2 ) = 




/pT|0i) + «1A 


/p~2\h), 




1^3) = 


«3V 


/pl|0 3 ) - «4A 


/Pl|04), 


(4-71) 


\xt) = 


— ia^ 




/Pl|04), 
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where 



(pi + P2) sin 20 + \j^p\p 2 + (p 4 - p 2 ) 2 sin 2 20 



oi 2 



\/2 ( Ap x p 2 cos 2 20 + (pi + p 2 ) 2 sin 2 26* + (p x + p 2 ) sin 26^Ap l p 2 + (pi - p 2 ) 2 sin 2 26* 

y/2p 1 p 2 cos 2$ 

4 Pl p 2 cos 2 20 + (pj + p 2 ) 2 sin 2 20 + (p x + p 2 ) sin 26^A Pl p 2 + (pi - p 2 ) 2 sin 2 20^ ' 

(4-72) 

(P3 + Pa) sin 26* + \jAp 3 p 4 + (p 3 - p 4 ) 2 sin 2 26 



a 3 = 



V2 4p 3 p 4 cos 2 26 + (p 3 + p 4 ) 2 sin 2 26 + (p 3 + p 4 ) sin 26J4p 3 p 4 + (p 3 - p 4 ) 2 sin 2 20 



\/2p 3 p 4 cos 20 

0:4 



4p 3 p 4 cos 2 20 + (p 3 + p 4 ) 2 sin 2 20 + (p 3 + p 4 ) sin 26^4p 3 p 4 + (p 3 - p 4 ) 2 sin 2 20 
Now it is easy to evaluate Aj which yields 

Ai = \ ((pi - p 2 ) sin 20 + ^4p 1 p 2 + (p 1 -p 2 ) 2 sin 2 20^ , 
A 2 = § ((p 2 - Pl ) sin 20 + y/A PlP2 + (p 1 -p 2 ) 2 sin 2 20^) , 
A 3 = I ((p 3 - P4) sin 20 + ^4p 3 p 4 + (p 3 -p 4 ) 2 sin 2 20^) , 
A 4 = § ( (p 4 - p 3 ) sin 20 + y/4p 3 p 4 + (p 3 - p 4 ) 2 sin 2 20 



(4-73) 



Therefore, coefficients Ki = f 1 are: 



K 4 



pi 

5 (ipi-P'l) sin 26+\J 4pip2 + (pi-p2) 2 sin 2 2#j 



K 2 - Pi 



K 



3 — ■ 

2 



I ((P2-pi) sin2e+ v /4pip 2 + (pi-P2) 2 sin 2 26 J' 

P3 

5 ((P3-P4) sin2e+ v /4^ 4 +(P3-P4) 2 sin 2 26»j ' 



(4-74) 



K* = —, 



\ {(Pi-P'i) sin 2e+-\/4p3P4+(p3-p4) 2 sin 2 

Writing the ICD state in wootters basis, we evaluate its robustness of entanglement with 
respect to the set separable state, diagonal in ICD basis, simply by chossing the separable 
states p' s and p" s on the corresponding boundaries, as follows, 

P's = \ ((Pi - V 2 ) sin 20 + v /4p i p 2 + (p / 1 -p 2 ) 2 sin 2 20) {x[ \x[) 



+ 2 
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\ ((P 2 -P'i) sin 20 + vW 2 + (pi-P2) 2 sin 2 20) 
+ ^ (Vs -P4) sin 20 + v/VsK+^s-K) 2 ^^^ 
+ \ ((P4 - Ps) sin 20 + yVsK + (Ps " P 4 ) 2 sin 2 2fl) (4|<) , 
P'i = ^ ((p? " P2) sin 20 + V / V/p 2 / + (p / i , -p 2 / ) 2 sin 2 20) (a^) (4-75) 
+5 ((P 2 " Pi) sin 20 + V /4KP , 2 ' + (P , 1 , -P , 2 , ) 2 sin 2 20) (x' 2 |x 2 > 
4 ((PS " P4) sin 29 + VWi + (P3 " Vl) 2 sin 2 20) ( : 



+- [ (//! - //,') sin 20 + y/ l//{//f + (,/<' - //[)-siir 20 j (JM 
+\ ((P'i - P3) sin 20 + v / 4p / K + (p , 3 , -P4) 2 sin 2 20) (x' 4 \x' A ) , 
with the corresponding robustness of entanglement as: 

mm(Ki + Kj) 



s = 



C (4-76) 



where are given in (4-74). 

It is obvious that, BD state correspond to particular case of = |, ICD ones with Ki — 1 
for % = 1,2,3,4. Therefore, the robustness of entanglement given in (4-76) become s = C 
which is agreement with (4-27) 

4.4 2 ® 3 Bell decomposable state 

In this subsection we obtain robustness of entanglement for Bell decomposable states of 2 <g> 3 
quantum systems. A Bell decomposable density matrix acting on 2 <g> 3 Hilbert space can be 
defined by 

P = X>I^>(^|, 0< Pi <l, J> = 1, ( 4 - 77 ) 

i=\ i=l 

where are Bell states in H & = H 2 ® if 3 Hilbert space, defined by: 



hM = ;^(|ll> + |22», |^> = -^(|11>-|22», 
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l^> = ^(112) + |23», m = -^(|12) - |23», (4-78) 
|^5> = ^(|13> + |21», |^ 6 ) = -^(|13)-|21». 

It is quite easy to see that the above states are orthogonal and hence it can span the Hilbert 
space of 2® 3 systems. From Peres- Horodeckis [34, 35] criterion for separability we deduce that 
the state given in Eq. (4-77) is separable if and only if the following inequalities are satisfied 

(Pi -P2? < (P3+P4)(P5 +Pe), (4-79) 

(p 3 -Pa) 2 < (P5+P6)(pi+P2), (4-80) 
(P5-Pe) 2 < (pi+P2)(P3+Pa)- (4-81) 

In the following we always assume without loss of generality that pi > p 2 , Pz > Pa and p 5 > p 6 . 
Now in order to obtain robustness of entanglement for BD state given in Eq. (4-77) we choose 
p' s = HiPi iV'i) (V'il an( i P = HiPi l^i) (^il- We also assume without loss of generality that p s 
lies on the separable-entangled boundary defined by (all other cases where p s lies on other 
surfaces can be treated similarly) 



Pi ~ P2 = vVs+^Ws + Pe)- (4-82) 

Moreover p s must satisfies the other two separability conditions (4-80) and (4-81). This means 
that entangled state p violates separability condition (4-79), i.e. we have 



P\>P2 J r\l (P3 + Pa) (P5 + Pa) ■ (4-83) 
Here the boundary of the separable states is given by 

((Pi -P2? - (Ps + P^iPn+Pe)) = 0. (4-84) 

States that saturate inequality (4- 79), (4-80), (4-81)) form a plane called S\ . All states 
violating these three inequality are entangled states. These states form an entangled region 
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with Si as its separable boundary. There exist, however, other entangled regions correspond 
to the dominating pj (j = 2, ..,6), respectively. These regions also define separable planes Sj. 
The Planes <Sj together with the planes <S ? ', correspond to p^ — 0. 

Below in the rest of this subsection we will use Eqs. (3-5) and (4-21) to calculate robustness 
of entanglement for 2 <g> 3 entangled Bell decomposable density matrix 

6 . 

P = ^2Pi\^i)(^i\ with pi > p 2 + \J(p3 + p 4 )(P5 +Ps), (4-85) 
i=i 

and an arbitrary separable density matrix 



Pa = S Pi I & I ' Pi ^ Pi + V (P3 + Pi) (P5 + Pe) ( 4 " 86 ) 
i=l 

The SDP optimization of —Tr(Ap) with respect to — Ap > yields 

A = mmAMA4% (4-87) 

Pi P2 P3 Pi PS P6 

The choice of A = ^ , i — 1, ...,6 consistent with positivity of p — sp^, implies that p" lies at 
the boundary <Sj , i = 1, 6. Numerical calculation indicates that the minimum s is achieved 
with the choice of A = ^ (for more details refer to ref[36]). Therefore the robustness of 
entanglement becomes 

Sl = (4-88) 

Pi 

we see that s is a monotonic decreasing function of p\ and its minimum is achieved for p\ = 



p' 2 + \J (p' 3 + p 4 ) (ps + pg) which implies that p' s lies at the boundary as separable state, too, in 
agreement with those results of ref[36]. 

Now in order to obtain minimum robustness of entanglement with respect to the set of 
separable states diagonal in 2 x 3, BD states, all we need to draw a line from p to interest the 
separable boundary 

Pl=P2 + yJip'a+p'M+p'e), (4-89) 

at p' s and the other boundary 



Pi = P2 + V(P3 + PZ) (P5 + Pe)- (4-90) 
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at p" , respectively. 

Therefore, from the relation p = (1 + s)p — sp" we have 

Pi = (1 + s)p\ - sp'i (4-91) 

Hence the robustness of entanglement s\ becomes 

s = (gi ~ P2) 2 ~ (P3 + P4)(p 5 + gg) / 4 g2 x 

1 1 -2(px - pa) (p? - p 2 ') + (p 3 + Pa) (pg + pg) + (p 5 + Pe) (pg + P4) ' 

where the maximization of the denominator, by using the Lagrange multipliers method due to 

existence of constrains (4-90) and normalization of Tr{p") = 1, leads to the following results 

for the minimums robustness of entanglement 

Si = 3((pi-p 2 ) 2 - (P3+P 4 )(P5+P6)) , 4 _ g3 j 

2(J(2 Pl - l) 2 + 3((px - p 2 ) 2 - (p 3 + p 4 ) (p 5 + Pe)) - (2p 2 - 1)) ' 



4.5 Werner states 

Werner states are the only states that are invariant under local unitary operations and for 
d <S> d systems the Werner states are defined by [37] 

where I stands for identity operator and T = ^ ^ s s b° wn that Werner state is 

separable iff < / < 1. 

Now to obtain the optimal robustness of entanglement with respect to accessible separable 
region, that is, density matrices of werner type with < /' < 1, all we need is to choose an 
arbitrary point p' s in separable region < /' < 1. Then the SDP method of optimization of 
—Tr(pA) with respect to p' s — Ap > yields 

A = mm {___ ) __ _} = __ (4-95) 

where the second equality follows from the fact that parameters / and /' are restricted to the 
regions / G ( — 1, 0) and /' G (0, 1), respectively. 
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Therefore, for the corresponding parameter s we get 

f - f 1 - f 

-1 + - 77, (4-96) 



-f + i i-r 

which is a monotonic decreasing function of /'. Hence, the optimal robustness of entanglement 
is s = —f which corresponds to the choice of werner type separability matrix with /' = 0. On 
the other hand we have 

which corresponds to the separable states with / = 1 in Eq.(4-94). 

Again both separable states p' s and p" are at the boundary of separable region in agreement 
with [1, 31, 3]. 

4.6 Isotropic states 

The d®d bipartite isotropic states are the only ones that are invariant under U®U* operations, 
where * denotes complex conjugation. The isotropic states of d<g>d systems are defined by [38] 

1 -F 



P 



d? 



— (i- 0<F<1, (4-98) 



where = ^ J2i N) is maximally entangled state. It is shown that isotropic state is 

separable when < F < | [38]. 

Again to obtain the optimal robustness of entanglement with respect to accessible separable 
region, that is, density matrices of Isotropic type with < F' < ^, all we need is to choose an 
arbitrary point p' s in separable region < F 1 < |. Then the SDP method of optimization of 
— Tr(pA) with respect to p' s — Ap > yields 

^ pi pi pi 

A = min{— -,—} = —, (4-99) 

where the second equality follows from the fact that parameters F and F' are restricted to the 
regions F e 1) and F' G (0, l/d), respectively. 
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Therefore, for the corresponding parameter s we get 

8 = y, ~ 1, (4-100) 

which is a monotonic decreasing function of F' . Hence, the optimal robustness of entanglement 
is s = dF — 1 which corresponds to the choice of Isotropic type separability matrix with F' — ^. 
On the other hand we have 

p" = p-(s + l)p' s = -^—^ (l - |^+) (^+|) , (4-101) 

which corresponds to separable states with F = in Eq.(4-98). 

Again both separable states p' s and p" are at the boundary of separable region in agreement 
with [1, 31, 3]. 



4.7 One parameter 3 3 state 

Finally let us consider a one parameter state acting on H 9 = H 3 <g> H 3 Hilbert space as [30] 

ct 5 — ol . 
+ + 2<a<5, (4-102) 

where 

l^ + ) = ^(|H) + |22) + |33)), 

a+ = | (|12) (12| |23) (23| + |31) (31|) , (4-103) 

<7_ = i(|21)<21||32) (32| + |13)(13|). 
p is separable iff 2 < a < 3, it is bound entangled iff 3 < a < 4 and it is distillable entangled 

state iff 4 < a < 5 [30]. 

Similarly in order to obtain the optimal robustness of entanglement with respect to acces- 
sible separable region, that is, density matrices of 3 <8> 3 type with 2 < a' < 3, all we need is 
to choose and arbitrary point p' s in separable region 2 < a' < 3. Then the SDP method of 
optimization of — Tr(pA) with respect to p' s — Ap > yields 

A = — , (4-104) 

a 
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where the second equality follows from the fact that parameters a and a' are restricted to the 
regions a G (3,5) and a' G (2,3), respectively. 

Therefore, for the corresponding parameter s we get 

s = - 1, (4-105) 

a 

which is a monotonic decreasing function of a'. Hence, the optimal robustness of entanglement 
is s = | — 1 which corresponds to the choice of 3 ® 3 type separability matrix with a' = 3. On 
the other hand we have 



2 



5 



p"=p-(s + l)p' s = - V + ) (^ + + (4-106) 

which corresponds to the separable states with a = in Eq. (4-102). 

Again both separable states p' s and p" are at the boundary of separable region in agreement 
with [1, 31, 3]. 

4.8 Multi partite isotropic states 

In this last subsection we obtain robustness of entanglement for a n-partite d-levels system. 
Let us consider the following mixture of completely random state p = 1 /d n and maximally 
entangled state |?/> + ) 

p = (1 - r )^- + r\ij + )(ij + \, 0<r<l, (4-107) 

where I denotes identity operator in <i n -dimensional Hilbert space and \^ + ) = Yh=i \ii • • - i)- 
The separability properties of the state (4-107) is considered in Ref. [39]. It is shown that the 
above state is separable iff0<r<r = (l + d n ~ v ) 1 . 

Again to obtain the optimal robustness of entanglement with respect to accessible separable 
region, that is, density matrices of multi partite isotropic type with < r' < r , all we need 
is to choose an arbitrary point p' s in separable region < r' < r . Then the SDP method of 
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optimization of — Tr(pA) with respect to p' s — Ap > yields 

A = i±lKLzi), (4 . 108) 

where the second equality follows from the fact that parameters r' and r are restricted to the 
regions r' G (0, r ) and r G (r , 1), respectively. 

Therefore, for the corresponding parameter s we get 

(r-r')(d n - 1) 



l + r'(d" - 1) ' 



(4-109) 



which is a monotonic decreasing function of r'. Hence, the optimal robustness of entanglement 
is s = ^+ro(i" J~i)^ which corresponds to the choice of multi partite isotropic type separability 
matrix with r' = r . On the other hand we have 

p" = p - (s + l)p' s = -^—^1 - |^ + ) (^ + |), (4-110) 

which corresponds to separable states with r = in Eq. (4-107). 

Again both separable state p' s and p" are at the boundary of separable region in agreement 
with [1, 31, 3]. 



5 Appendix 

In the definition of the robustness we have to minimize over all separable states. But almost in 
all of examples, we have considered a specific set of separable states, i.e, the diagonal separable 
states (in the basis that the entangled state itself is diagonal). Here in this appendix we try to 
show that the minimum robustness of an entangled diagonal density matrix in a given basis, 
with respect to the set of separable diagonal states in the same basis, is also minimum over 
the separable sets of off-diagonal extension of these diagonal separable states. In [1], it has 
been shown that the robustness given in equation (4-17) is minimum over all Bell decomposable 
states. In the following, we generalize it and show that for any orthogonal basis the off-diagonal 
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elements of p' s and p" (in basis that p is diagonal) play no role in robustness. First we consider 
the diagonal entangled density matrices in some orthonormal basis and we consider the generic 
separable states p' s and p" defined by 



P's = J2p'i\Xi)(Xi\ +J2 a iMi)(Xj\ (5-1H) 

i i,j 

p" = Y,Pi\Xi)(Xi\+Y,h\Xi)(Xj\, (5-H2) 
where an = bu = and \xi), i = 1, 2, .. are orthogonal states. Then the pseudomixture equation 



(4-23) implies that the following equations must hold 

Pi = (1 + s)p[ - sp", 



[1 + s)a,ij — sbij = 0. 



Now one can easily obtain robustness of p relative to p' s as 

" P - P's II 



(5-113) 
(5-114) 



s = 



P's - P" 



E,fa-P0 2 + Tr(^) 
\EM-p'{) 2 + Tr((A-By) 



E,fa-P0 2 + Tr(^) 

where A and B are Hermitian matrices with non vanishing off-diagonal matrix elements 
and bij, respectively, and in the last line we have used equation (5-114). By solving equation 
(5-115) for robustness s we get 

'^fc-ri) 2 (5-116) 



Equation (5-115) shows that off-diagonal elements of p' s and p" (in basis that p is diagonal ) 
play no role in robustness. Now we consider the general case of diagonal entangled density 
matrices in some non-orthogonal basis. Obviously in this case the entangled density matrix p 
can be written as 

P = ^2Pih)(Vi\ = ^2Pij\f}i)(i}i\, (5-H7) 
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where \r)i),i = 1,2, ... are non-orthogonal states and \fji),i = 1,2, .. are dual of \r)i),i = 1,2, ... 
, i.e., we have (fji\\r)j) = Sij,i,j = 1,2, .. . Also the separable states p' s and p" can be written 

as 



\Vi)(Vi\ +J2 a V 




= 

i 


Vi)(Vi\ + E fi *i 




(5-118) 


\Vi)(Vi\ + E% 


Vi)(Vj\ 


i 


\Vi)(Vi\ 


Vi)(Vj\, 


(5-119) 



where an = bu = an = bn = 0. again the pseudomixture equation (4-23) implies that the 
following equations must hold 

Pi = (l + s)p' t - sp'i (5-120) 

= (1 + 8)0*3 ~ sh iv (5-121) 

pa = (1 + s)p' t - sp" t , (5-122) 

Pij = (1 + s)a,ij - sbij. (5-123) 
Now the robustness of p relative to p' s can be easily obtained by using above Equations as 



P~ P's 



P's - P" 



Ei{Pi-p'i)(Pu-p'i) / 5 124 n 

Again this equation shows that off-diagonal elements of p' s and p" in non-orthogonal basis play 
no role in robustness. 



6 conclusion 

Using the elegant method of convex semidefinite optimization method, we have been able to 
obtain the robustness of some set of mixed density matrices with respect to some accessible 
separable set. In this method we have been able to calculate the robustness without using any 
kind of the space of density matrices, where the results that obtained are in agreement with 
those of norm-method of ref[2]. 
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Also using SDP method we have shown that the separable density matrices contributing 
to the robustness lie at the boundary of accessible separable region. 
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Figure 1: All BD states are defined as points interior to tetrahedral. Vertices Pi, P 2 , P3 and 
P 4 denote projectors corresponding to Bell states given in Eqs. (4-8) to (4-11), respectively. 
Octahedral corresponds to separable states. 

Figure 2: Entangled tetrahedral corresponding to singlet state. Point t denotes a generic 
state p and points t' and t" are on the separable boundary planes. 

Figure 3: The space of a generic two qubit density matrix is represented by a tetrahedral. 
Vertices P for % = 1,2,3,4 correspond to pure states defined by p = A; Irregular 
octahedral corresponds to separable states. Separable planes Si and S[ are shown explicitly. 



